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Abstract 

We give 50 digits values of the simple continued fractions whose denomi- 
nators are formed from a) prime numbers, b) twin primes, c) generalized 
d-twins, d) primes of the form -|- n'^, e)primes of the form -|- 1, 
f) Mersenne primes and g) primorial primes. All these continued frac- 
tions belong to the set of measure zero of exceptions to the theorems 
of Khinchin and Levy. We claim that all these continued fractions are 
transcendental numbers. Next we propose the conjecture which indi- 
cates the way to deduce the transcendence of some continued fractions 
from transcendence of another ones. 



1 Introduction 

Let oo be an integer and let Ofc, k = 1,2, . . . ,n are positive integers (in general 
can be arbitrary complex numbers, see e.g. [30]). Then 



r = [ao; ai, 02,03, ... ,a„] = ao H (1) 
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is the simple (i.e. with all nominators equal to 1) finite continued fraction. The 
numbers a^, k = 1,2, . . . ,n are called partial quotients and 

— = [ao;ai,a2,a3, . . . ,ak], k = l,2,...n (2) 

is called the k-th convergent of r. If for the infinite continued fraction 

[ao; 01,02,03,...] (3) 

the sequence of convergents Pn/Qn converges to some limit r when n — )■ oo then 
we say that the infinite continued fraction [oq; Oi, 02, 03, . . .] is equal to r. The 
convergence of the continued fraction ([s]) is linked to the behavior of the sum of 
partial quotients o„: 

p °° 
sequence — ^ is convergent to r o„ is divergent (4) 

^" n=l 

see e.g. [291 Theorem 10, p. 10]. It means that for convergence of the continued 
fraction it is necessary that both P„, Qn — ?■ 00 in such a way, that the ratio Pn/Qn 
has a definite limit for — )■ 00. If the infinite continued fraction is convergent then 
the values of the convergents -Pfc(r)/(5fc(r) approximate the value of r with accuracy 
at least 1/QkQk+i [29l Theorem 9, p. 9]: 



Qk 



< ^ < ^ < J_ (5) 

QkQk+l Q\cLk+l Qk 



Rational numbers have finite continued fractions, quadratic irrationals have periodic 
infinite continued fractions and vice versa: eventually periodic continued fractions 
represent quadratic surds. All remaining irrational numbers have non-periodic con- 
tinued fractions. 

Khinchin has proved that [29^ p. 93] 

log2 m 



lim (01O2 . . . On) " = n S 1 + ^ > =Ko^ 2.685452001 . . . 

n-s>oo ^ ^ 1^ m(m + 2) J 



(6) 



is a constant for almost all real r, see also [39], [231 §l-8]- The exceptions are rational 
numbers, quadratic irrationals and some irrational numbers too, like for example the 



Euler constant e = 2.7182818285. . ., but this set of exceptions is of the Lebesgue 
measure zero. The constant Kq is called the Khinchin constant. 

In 1935 Khinchin [28] has proved that for almost all real r the denomina- 
tors Qn{r) of the convergents of the continued fraction expansions for r satisfy 
(r) = Lq and in 1936 Paul Levy [32] found an exphcit expression for 
this constant Lq: 

lim ^Qn{r) = e"'/^2iog(2) ^ = 3.27582291872 ... (7) 

n— >oo 

All presented below continued fractions belong to this exceptional set of irrationals 
for which the geometric means of the denominators {aia2 ■ ■ ■ CLnY^'^ and the n-th 

1 /n 

radical roots of the denominator Qn tend to infinity. 

In this paper we will consider continued fractions with partial quotients given by 
an infinity of all primes as well as primes of special form belonging to families contain- 
ing conjectured infinity of members. All these continued fractions are non-periodic, 
and thus are irrational, but we claim that all of them are also transcendental. In 
Sect. 3 we review some facts and theorems concerning the transcendentality of con- 
tinued fractions. In Sect. 4 some experimental results regarding transcendentality 
of numbers constructed from primes are presented. 

2 Seven examples 

In consecutive sections we will discuss the following cases: the set of all primes 
2, 3, 5, 7, . . ., twin primes, generahzed d-twins, i.e. pairs of adjacent primes separated 
by d, primes of the form form + n"^, primes given by the quadratic form + 1. 
Next are considered sparse Mersenne primes and at the end even sparser primorials 
primes. 

It is possible to consider other families of primes, like Sophie Germain primes (it 
is conjectured that there are infinitely many of them), irregular primes (Jensen in 
1915 proved that there are infinitely many of them), regular primes of which it was 
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conjectured that e^^^^ ~ 61% of all prime numbers are regular, the CuUen numbers 
n2^ + 1 when they are primes etc. but we leave it for further studies. 

Except Sections 2.1 and 2.4, where we will treat the families of primes containing 
rigorously proved infinity of members, all remaining consideration are performed 
under the assumption there is infinity of primes in each class of primes, although 
proofs of infinitude of all these sets of primes seems to be very far away. Thus many 
of our reasonings are heuristical. 

The examples are in order of sparseness of each family of primes. 

2.1 The set of all primes 

Let us put a„ = p„ where p„ denotes the n-th primes: [0; 2, 3, 5, 7, 11, 13, . . .]. As 
there is an infinity of primes the condition Q is fulfilled and let us denote the limit 
of the continued fraction by 

1 

u = [0; 2, 3, 5, 7, 11, 13, . . .] = (8) 



1 



2 + 



3 + 



5+ 

1 

7+ 

11+ ■•• 

Using PARI system and all 1229 primes up to 10000 it is possible to obtain over 

8000 digits of the above continued fraction in just a few seconds because 

, 3.38592889 . . . x 10^^97 
|0;2,3,5,7,ll,13,...,99731 = ^^3^^^^5^-^ (9) 

and the product of QkQk+i on the rhs of ([s]) is larger than 10^^°°. The first 50 digits 
of u reads: 

u = 0.43233208718590286890925379324199996370511089688 .... (10) 

This number is not recognized at the Symbolic Inverse Calculator (http:/ /pi.lacim.uqam.ca/eng/) 
maintained by Simone Plouffe. Accidentaly, it is very close to the one of Renyi's 



parking constants = (1 - e-'^)/2 = 0.43233235838..., see ^ pp. 278-283]: 
mR-M = 2.712... X 10~^ 

It is possible to obtain analytically the geometrical means of the partial quotients 
in ([s]). It is well known (see e.g. [211 Chap. 4]), that the Chebyshev function 9(x) 
behaves like: 

e{x)=J2^0g{p)=X + 0{^/^). (11) 

Thus skipping the error term we have 



llpk = eP-. (12) 



k=l 



It is well known that ^37t Sect. 2. II. A] that 



Tl log log(7?/) 

Pn = nlog(n) + n(loglog(n) — 1) + o ( ) . (13) 



n 



For our purposes it suffices to know that 

Pn > n\og{n) for n > 1 (14) 

see e.g. [38]. Hence we can write for the geometrical means of the partial quotients 
the estimation: 

(n \ " 

Y[pk\ ={eP")^>n^oo (15) 

thus the continued fraction u belongs to the set of measure zero of exceptions to the 
Khinchin Theorem ([6]). It is also an exception to the Levy Theorem, because from 
the general properties of continued fractions: 

Qn+l = O'nQn + Qn-l (16) 

we have Qn > 11^=1 Pfc > and thus QI/"^ — )■ oo in contrast to (7). It is an explicit 
example of the continued fraction with unbounded {QnY^"'- 
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2.2 Twin primes 

The twin prime conjecture states that there are infinitely many pairs of primes 
differing by two: — t„ = 2. Let 712(0;) denote the number of pairs 

of twin primes {tn,tn+i) smaller than x. Then the conjecture B of Hardy and 
Littlewood [26] on the number of prime pairs p,p + d applied to the case d = 2 gives, 
that 

u X 

7C2{X)^C2 2i ^ du = C2 ^ +..., (17) 

J2 log [u) log [x) 

where C2 is called "twin constant" and is defined by the following infinite product: 

^2 = 217(1-- — ^1= 1.32032363169... (18) 

If there is indeed (as everybody believes, see e.g. an infinity of twins, then the 
continued fraction 

U2 = [0; 3, 5, 5, 7, 11, 13, 17, 19,...] (19) 

should be infinite, non-periodic and convergent. We count here 5 two times as it is 
a customary way of defining the Brun's constant |13] and it an only case of double 
appearance of a prime in the set of twins as for adjacent twin pairs {p — 2,p) and 
{p, p + 2) one of numbers {p — 2,p,p + 2) always is divisible by 3. Again performing 
calculations in PARI and using primes < 10000 we found here 205 twin pairs (but 



only 409 different primes) and first 50 digits of the continued fraction (19) are 



U2 = 0.31323308098694591263078648647217280043925117451 .... (20) 

There is much less terms in U2 up to 10000 than primes < 10000 in u, hence the value 
of U2 was obtained with accuracy about 2900 digits. We have checked using Plouffe's 
Symbolic Inverse Calculator (http://pi.lacim.uqam.ca/eng/), that this constant is 
not recognized as a combination of other mathematical quantities. 

Because twin primes are sparser than all primes we have tn > Pn thus in view 



of (15) the geometrical means (3 5. . .t„)i/'^ will diverge even faster, hence the 



continued fraction U2 belongs to the set of exceptions to the Khinchin Theorem. It 



is also a counterexample to the Levy Theorem, because denominator Qn{u2) of the 
n-th convergent of U2 is larger than the denominator Qn{u) of the n-th convergent 
of u. 

2.3 Generalized (i- twins 

It is natural to consider the whole family of continued fractions Ud-, c? = 2, 4, 6, 8, . . . 
formed from the consecutive primes separated by d: Pn+i — Pn = d. We put this 
example here after twins, although for sufficiently large d the primes Pn+i — Pn = d 
will be even sparser than say Mersenne primes and from the other side d = 2 are 
less frequent than d = 6, see [35j. The consecutive primes separated by d = 4 are 
sometimes called Cousins, [50]. For example, in the case of d = 6 we have : 

U6 = [0; 23, 29, 31, 37, 47, 53, 53, 59, 61, 67, 73, 79, 83, 89, 131, 137, 151, 157, 157, 163, ...] 

and some primes p„ when p„ — Pn-i = Pn+i — Pn do appear twice (in the case of 
U2 only 5 appears two times). As in the case of twins it is conjectured that for 
each d there is an infinity of prime pairs {pn+i^Pn) with Pn+i — Pn = d, see e.g. 
[To] . [35j . From this conjecture it follows that the numbers Ud are irrational. Using 
PARI/GP we have calculated the values of Ud up to d = 570, what took four days 
of CPU time on the 64 bits AMD Opteron 2700 MHz processor. We have searched 
for primes up to 2^^ = 1.759. . . x 10^'^ and the largest encountered gap between 
consecutive primes was d = 706 which appeared only once. We have calculated Ud 
if there was a number of gaps of given d sufficient to determine Ud with at least a 
few hundreds digits (for example, there were 17 pairs of 570-twins up to 2^^). The 
Table I gives a sample of obtained values with 50 digits accuracy. The whole file 
with 275 values of Ud given with more than 110 digits is available from the author 
webpage http:/ /www. ift.uni.wroc.pl/~mwolf/u_d.dat. 

For large d the value of Ud is practically determined by the first occurrence Pf{d) 
of that gap — pairs of consecutive primes with gap d ^ 2 are separated by very large 
intervals, for example first d = 540 appears between (738832927927, 738832928467) 
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and next gap d = 540 is between (3674657545087,3674657545627). It was conjec- 
tured by Shanks in 1964 [42] that the gap d appears for the first time at the prime 
Pf{d) ~ e^. We have given heuristic arguments [IHj that 

Pf{d) ~ v^exp Q0n2(rf) +4d^ (21) 

and for ^ 1 simply Pf{d) ~ \fde'^. Thus we claim that for large d there should 
be the approximate formula: 

[0;v^e^,v^e^ + c/] ^ — i^. (22) 

The plot of and comparison with the Shanks and our conjecture is given in 
the Fig.l. 

Again like u and continued fractions Ud belongs to the set of exceptions to 
the Khinchin Theorem and Levy Theorem. 

2.4 Primes of the form rr? + 

In the seminal paper [24J John Friedlander and Henryk Iwaniec have proved that 
there exists infinity of primes of the form -|-n^. More precisely, if 7rFi(x) denotes 
the number of primes of the form + < x then approximately 

vrFi(x) ~ (23) 
log(x) 

where the constant C^i = V2T{lf /3tt^/^ = 1.112835788898764 ... and here T is the 
Euler Gamma function. Thus taking as partial quotients of the continued fraction 
primes of the form + n'^ for sure we will obtain an irrational number which we 
will denote upi. 

UFi = [0; 2, 5, 17, 17, 37, 41, 97, 97, . . .] (24) 

Like in previous examples some primes appear twice: 17 = 4^ + I"' = 1^ + 2^, 97 = 
9^ + 2^ = 4^ + 3'^ etc. Looking for all primes of this form with 1 < m < 100 and 
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Fig.l The plot of Ud and two approximations: in green the Shank's conjecture 
1 1 



and in red our conjecture 



1 



1 



Vde-^ + d 



1 < < 10 (the largest prime was 19801 = 99^ + 10^) we get the value of upi with 
over 1100 digits accuracy; the first 50 digits of it are: 

UFi = 0.455024816490170022369052808279744824105755548905 . . . (25) 

Let us notice that 

1/(2 + 1/(5 + 1/(17 + 1/(17 + 1/(37 + 1/(41 + 1/(97 + 1/98))))))) = 
20993638525 

^—-^ = 0.455024816490170022369048157801049432084768331968 . . . 

46137348479 

and the difference between this value and upi is less than 10"^'^ ! 
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Table I 



d 


Ud 


4 


1.4103814184127409729946079947661391024642878552250x10-^ 


6 


4.3413245800886640441937906138426444157119875018764x10-2 


8 


1.1234653732060451418609230935360294984983811524705x10-2 


10 


7.1938972705064358418419102215951120335820544247877x10-3 


12 


5.0250059564863844924667112008186998625931272954692x10-3 


14 


8.8489409307271044901495673780577102976304420791245 x 10-^ 


16 


5.4614948350881467294308534284337241698766002935218x10-^ 


18 


1.9120391314299159400657740968697274305281924125799x10-3 


20 


1.1273943145526585257207207582991176443515379616999x10-3 


22 


8.8573891094929851372874303656530678911673854053699x10-^ 


24 


5.9916096230989554005997263265407846890656053212565x10-^ 


26 


4.0371410525148524468010569219212401713453876041188x10-^ 


28 


3.3658696996531260967017397551173798914121748535404x10-^ 


30 


2.3272049015980164345521674554989676374011829679698x10-^ 


32 


1.7885887465418665415382499015390795012182483537844x10-^ 


34 


7.5357908538425634007656299916322144807040843935028 x 10"^ 


36 


1.0470107727765143055064789951193220804598138780293x10-^ 


38 


3.2687215994929278130910770751451289367042590019431 x 10-^ 


40 


5.1725029603788623137563671868924142637218125718293x10-^ 


42 


6.1954029872477528100134249220879079074519481392595x10-^ 


44 


6.3763310332564890009355447509046689625278954819441 x lO'^ 


46 


1.2275511580096446939755547564625149207372813752259x10-^ 


48 


3.5424563347877377649245656903453981296411399487963x10-^ 


50 


3.1341084997626641267187094247975857118584579732840x10-^ 






566 


2.0417988154535953561248601983565125430801657124094x10-13 


568 


1.6638019955234637865242752590874891539355008513604x10-13 


570 


2.2511824714719308536000694530283450847909292429681 x 10-^3 
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2.5 Primes of the form + 1 

Now let us consider the set of prime numbers 



Q = {2, 5, 17, 37, 101, 197, 257, 401, 577, 677, 1297, 1601, . . .} (26) 

given by the quadratic polynomial + 1 and let g„ denote the n-th prime of this 
form. By the conjecture E of Hardy and Littlewood [2S] the number 7Tq{x) of primes 
Qn < X of the form g„ = + 1 is given by 



log(x 

where 



^.(^) ~ C,-^, (27) 



/ (_l)(p-l)/2\ 

= 11 M - ^ '— — j = 1.372813462818246009112192696727 ... (2^ 



p>3 

.2 



Comparing it with (23) we see that indeed primes m + 1 are sparser than primes 



+ n^. For example up to 10*^ there are 65162 primes of the form + and 



only 841 primes of the form m + 1. Although the conjecture (27) remains unproved 
there is no doubt in its validity. Thus let us create the presumedly infinite continued 
fraction by identifying a„ = > 1: 

Ug = [0; 2, 5, 17, 37, 101, 197, 257, 401, 577, 677, 1297, 1601, . . .]. (29) 

Using 841 primes of the form m^ + 1 smaller than 10* and performing the calculations 
in PARI with precision set to 20000 digits we get over 11000 digits of Ug as the ratio 
on the rhs of d5| was < 10^^^™'^. First 50 digits of Ug reads: 



Ug = 0.45502569980199468718020210263808421898137687948 . . . (30) 



Let us remark that upi — Ug = 8.833 ... x 10 



There is no known formula analogous to (11 ) for primes of the form ra + 1, but 



because g„ > p„ the geometrical means of 2 ■ 5 ■ 17 . . . g^i will diverge faster than ( 15 ) 
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It is possible to obtain very rough speed of divergence of (2-5-17. . .g„)^/". Namely, 



making use of (27) and inverting 7rg(g„) = n we get: 



2nlog(n/C„)V r^, f ''^ \ , i f''^ \ /^-,\ 
^ ^ +21og — log log — (31) 



Cq J \ Cq J \ Cq 



Because 2 > Cg it follows that 2 • 5 ■ 17 . . . g„ grows faster than 2^"(n!)^/C|" > (n!)^ 
and the Stirling formula for n! gives that (2 ■ 5 ■ 17 . . . QnY^^ grows faster than and 
again Uq is the exception to the Khinchin Theorem as well as to the Levy Theorem. 



2.6 Mersenne primes 

The Mersenne primes are the primes of the form 2^ — 1 where p must be a 
prime, see e.g. [371 Sect. 2. VII]. Only 47 primes of this form are currently known, 
see Great Internet Mersenne Prime Search (GIMPS) at www.mersenne.org. For 
many years the largest known primes are the Mersenne primes, as the Lucas-Lehmer 
primality test (applicable only to Ain = 2^ — 1) needs just a multiple of p steps, 
thus the complexity of checking primality of Ain is 0{log{Ain))- Let us remark that 
algorithm of Agrawal, Kayal and Saxena (AKS) for general prime p works in about 
0{log^'^{p)) steps and modification by Lenstra and Pomerance in about 0(log^(p)) 
steps. 

Again there is no proof of the infinitude of but a common belief is that 
as there are presumedly infinitely many even perfect numbers thus there is also an 
infinity of Mersenne primes. 

Let us define the supposedly infinite and convergent continued fraction um by 
taking a„ = A^„: 

um = [0;3, 7, 31, 127, 8191, 131071, 524287, 2147483647, . . .] (32) 

Using all 47 Mersenne primes 3, 7, 31, . . . , 2*^^^^^^°^ — 1 in a couple of minutes we 
have calculated um with the precision better than 10^121949117. ^^.^^ fjigj^^g Qf 
are: 

Um = 0.31824815840584486942596202748140694243806236564 . . . (33) 

12 
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16 



12 - 



8- 



o o o log log(M„) 
«e-nog(2)-log log(2) 

Fit: 0.385 n + 0.6691 




50 



Fig. 2 The plot of loglog(7Wn) and the Wagstaff conjecture (34). The fit was made to all 
known Mn and it is 0.3854n + 0.6691, while ne"'^ log(2)n - log log(2) 0.3892n + 0.3665. 
The rather good coincidence of loglog(A^n) and ( p4| ) is seeming, as to get original A^n's 
the errors are amplified to huge values by double exponentiation. 



Of course um is also the exception to the Khinchin and Levy Theorems in view 
of the very fast growth of um — Wagstaff conjectured [l8],that M.^ grow doubly 
exponentially: 

log2log2A^„ ~ ne"^. (34) 

where 7 = 0.57721566... is the Euler-Mascheroni constant. In the Fig. 2 we 
compare the Wagstaff conjecture with all 47 presently known Mersenne primes. 

2.7 Primorial primes 

If p„ is the 72— th prime number then numbers of the form 2x3x5---xp„= p^jj are 
called primorials and % stands here by analogy of exclamation mark in the factorial. 
The primorials are expressed directly by the Chebyshev function Q[x): 

p^^ = e^(P") = e(i+°(i))P'\ (35) 
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For some primes r„ the numbers r^jj ± 1 are primes. They are called primorial 
primes and are even sparser than Mersenne primes as we will see below. Despite 
this rareness of primorial primes it was conjectured that there is infinity of them 
[TB] and thus the continued fractions 

= [0;3,7,31,211,...] (36) 

obtained from primorial primes of the form r^^ + l will be at least irrational number, 
as well as the continued fraction build from primorial primes of the form r^jj — 1: 

= [0;5,29,2309,...]. (37) 

From the known presently only 22 (see sequence A005234 in OEIS) primorial primes 
r„tl + 1 we get the continued fraction 

Ur+ = [0; 3, 7, 31, 211, . . . , 42209^ + 1, 14582311 + 1, 36643911 + 1, 39211311 + 1] = (38) 

0.318248165083690124777685589996787844788657122331533049467 . . . 

with the error less than IQ"^^'^^''*^. Let us remark that Ur+ — um = 6.678 . . . x 10^^, 
although only three first primes (3, 7, 37) are the same. 

From all 18 presently known (see sequence A006794 in OEIS) primes of the form 
r„ll - 1 we get 

Ur- = [0; 5, 29, 2309, . . . , 1587711 - 1] = (39) 

0.198630157303503810875201233614346862875870630898479777625647 . . . 

with the error less than 10"^^'^^^. Chris K. Caldwell and Yves Gallot gave heuris- 
tic arguments [13j that there is infinity of primorial primes of both kinds. More 
precisely, they claim that the expected numbers of primorial primes of each of the 
forms rH ± 1 with r < x are both approximately e'''log(x). From n = e'^log(r„) we 
get that 

rn~e"/^\ (40) 
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where r„ stands for n-th prime giving the primorial prime of the form r^jj ± 1. Then 
the n-th primorial prime, and hence of Ur±, will be 

an = = exp I e^'^ ^ p I . (41) 

From the formula 

E/(P)= r /^,dx+f{2)U{2)+f{n){7r{n)-U{n))- r{7r{x)-U{x)}dx (42) 

p<n ^Og[X) J2 

where 7i{x) is the number of primes < x and li{x) = dx/ log(a;) is the logarithmic 
integral, see [3H1 eq.(2.26)] we get 

p = liin^) + error = — — -— + error'. (43) 
^ ^ ^ 21ogn ^ ^ 

p<n '^^ ' 



Here error' besides expressions on rhs in (42) contains also higher terms coming 
from the asymptotic expansion of li{x^\ Finally we obtain 



~ exp exp (2e ^n)^ 



(44) 



From this it follows that the n-th primorial prime is much larger than the n-th 
Mersenne prime ~ 2?^^" . Indeed, the ratio: 

log(r„tt) _ 1 /e^y ^ ^^^^ 



log(7W„) 21og(2)nV2 
grows with n. 



3 Continued fractions and transcendence 

There is a vast literature concerning the transcendentality of continued fractions. 
The Theorem of H. Davenport and K.F. Roth [TB] asserts, that if the denominators 
Qn of convergents of the continued fraction r = [ao; ai, a2, . . .] fulfill 

^log(n) log(log(Q„(r))) 

lim sup = oo (46) 

n n 

15 



then r is transcendental. This theorem requires for the transcendence of r very fast 
increase of denominators of the convergents: at least doubly exponential growth is 



required for |46| The set of continued fractions which can satisfy the Theorem of 
H. Davenport and K.F. Roth is of measure zero, as it follows from the Theorem 31 
from the Khinchin's book [29], which asserts there exists an absolute constant B 
such that for almost all real numbers r and sufficiently large n the denominators of 
its continued fractions satisfy: 

Qn{r) < (47) 

The paper of A. Baker [6] from 1962 contains a few theorems on the transcen- 
dentality of Maillet type continued fractions [31], i. e. continued fractions with 



bounded partial quotients which have transcencendental values. In the paper [3j B. 



Adamczewski and Y. Bugeaud, among others, have improved (46) to the form: 



Tf 1- log(log(Q„(r))) 

„'^Pn2/3(log(Q„(r)))2/3 1og(log(Q„(r))) 

then r is transcendental. 

Besides Maillet continued fractions there are some specific families of other con- 
tinued fractions of which it is known that they are transcendental. In the papers 
[HS] , [1| it was proved that the Thue-Morse continued fractions with bounded par- 
tial quotients are transcendental. Quite recently there appeared the preprint [2] 
where the transcendence of the Rosen continued fractions was established. For more 
examples see [5]. 

Taking as the partial quotients a„ different sequences of numbers leads to real 
numbers which very often turn out be transcendental. For example, the continued 
fraction s for which a„ = n: 

s= [0; 1,2,3,4,...] = 1/(1 + 1/(2 + 1/(3 + 1/(4 + ...)))) = 0.697774657964... (49) 

is transcendental. Let us mention that the continued fraction s' with all partial 
quotients equal to consecutive odd numbers: 

s' = [0; 1, 3, 5, 7, 9, . . .] = -i tan(i) = ^ ~ = 0.761594155955764888119 . . . (50) 
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is also transcendental. All these facts are special cases of the results obtained by Carl 
Ludwig Siegel in 1929 in a long paper [H]. In particular, the continued fractions 
sd = [0; cti, (32, ct3, ■ • ■] are transcendental when the a„'s are rational and form an 
arithmetical sequence of the difference D and first element A : = A + kD, k = 



1, 2, 3, . . .. Siegel mentioned exphcitly the continued fraction (49), see [15 or p. 231 
in Gesammelte Ahhandlungen vol. I]. He obtained these results as corollaries from 
the continued fraction expansion of the ratio of Bessel's functions J\{x) (see also [H 
formula 9.1.73]): 

/yf'f = + (51) 

Ja(2zx) X . ^ ' 

A + 1 , ^ 



X 

A + 2 



+ 



X A + 3 _^ •. 

X 

which Siegel has shown to be transcendental for rational A and algebraic x and 
where the Bessel function of first order is given by 

^^W = E .,„,r(jA + i) (2) • <^2' 

where r(x) is the Gamma function, see jH', first formula on p. 231 in Gesammelte 
Ahhandlungen vol. I] or [1, formula 9.1.10]. For A = and x = 1 and taking into 
account the relation J_„(x) = (— l)"J„(x), see e.g. [H formula 9.1.5], we get the 



value of the continued fraction (49): 



s=[0;l,2,3,4,...] = -2-^. (53) 



The awkward form (51 ) can be written in more pleasant form in terms of modified 



Bessel functions of the first kind lv[x) defined by the series: 



There is a following relation between l^ix) and Ju{ix): 

h{x) = {-tyj,{tx), (55) 
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2k 



see [U formula 9.6.3]. Writing A = X/x, D = 1/x, i.e. A = A/D, x = 1/D we turn 

Ia/d-1 (I) 



(51) to the more elegant form 

[A;A + D,A + 2D,...,A + nD, 



lA/D 



(56) 



For A = —1, D = 2 (oi for A = — 1/2 and a; = 1/2 in (51)) we obtain the value 
of the continued fraction s' defined by the formula (ISOj): 



/-3/2(i: 



s' = [0;l,3,5,7,9,...] = l + — ^ = 1+2— ^ 

-'-1/2(1; j_i 



3/2 ( 



The transcendence of (50) follows for x 



already to Lambert and Euler |22j : 



tan(a;) 



/2 
X 



(57) 



1) from the formula known 



(5^ 



X 



1 - 



X 



3- 



and the fact that tan(x) takes transcendental values at algebraic arguments. From 



this and from (57) as a byproduct we have the identity: 



tan(i) 



3/2 1 



(59) 



Another possibility for partial quotients is the geometrical series : an = q"" and 
we believe that the corresponding continued fractions: 



(60) 



are transcendental for all natural q > 2. This continued fraction is linked to the 
famous Rogers-Ramanujan continued fraction defined by 

.,1/5 



RRiq) 



(61) 



1 + 



q 



1 + 



1 + 
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From the general transformations of continued fractions rules see (SDJ P-9] we have 
the relation: 

q 1 

q'/'^RRiq) = 



1 + q-^ + 



q' 



q-^ + 



1 + + 



1 + q~' + 

q^ 1 
1 + — q-^ + 



+ ••. 

(62) 



In [19] [20] it was proved that RR{q) is transcendental for all algebraic |g| < 1, 
but it needs some further work to infer from this the transcendence of Gq. Let us 
mention, that quite recently K. Dilcher and K. B. Stolarsky [18] have proved that 
the continued fractions 

^9 = (63) 



are transcendental for all integer q > 2 — it follows immediately from (46) and 
the double exponential growth: Qnidq) > 11^=1 Q'^'^- Adamczewski [2j extended 
this to all complex |g| > 1 which are algebraic numbers. Another (family) class of 
transcendental continued fractions can be found in [T7j. 

Next we can construct a number / where partial quotients are factorials a„ = n!: 

/ = [0; 1, 2, 6, 24, 120, ...,n\,...] = 0.6840959001066225003396337 . . . (64) 

Even these partial quotients increase too slowly to apply the Theorem of Adam- 



czewski and Bugeaud (48). For large n we have approximately Qn{f) ~ n!fc=i^'- 
This last product is called superfactorial and denoted by n$, see also [25j, exercise 
4.55]. We prefer the notation n\' = YYk=i Superfactorial can be expressed by the 
Barnes G-function for complex z defined by 



G{z + 1) = (2vr)^/2e-(^(^+^)+^^')/2 J] [(l + ^)%--+-V2„j _ ^g^^ 



n=l 
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It satisfies the functional equation 



G{z + 1) = T{z)G{z) (66) 

and from this we have that 

n\-=G{n-2). (67) 
The analog of the Stirling formula for G{z) gives |17] : 

log G{z + l) = z'( - ? V J log(2vr) - • . • (68) 



AJ 2 
From this we obtain 

and unfortunately 

v/log(Q„(/))loglog(g.(/)) ^ ^ 

hence we do not get transcendentality of / via the Theorem of Adamczewski and 
Bugeaud. 

The continued fraction build from Fibonacci numbers a„ = Fn 

F = [0; 1, 1, 2, 3, 5, 8, ... , F„, . . .] = 0.588873952548933507671231121246787384 .... 

(71) 

appears at the Sloane The On-Line Encyclopedia of Integer Sequences as the entry 
A073822. 

Apparently both / and F also should be transcendental, but we are not aware 
of the proof of this fact. The factorial over Fibonacci numbers behaves as 

n 

Y[Fk = 0"{"+l)/25-n/2c< + O (0"{'^-3)/25-n/2) (72) 
k=l 

where = (1 + ^/5)/2 and G ^ 1.226742, see [2S1 Exercise 9.41] and it is too slow 
to use the Davenport - Roth Theorem. 

Let us quote at the end of this Section the following remarks from the |7l p. 
104]: "And the latter recalls to mind another outstanding question in Diophantine 
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approximation, namely whether every continued fraction with unbounded partial 
quotients is necessarily transcendental; this too seems very difficult" . Now there is 
a common believe that also algebraic numbers of degree > 3 have unbounded partial 
quotients, see e.g. [H], [3]. 

4 Transcendence of u, U2, Ud, u-pi, Ug, um, Ur± 

Because all considered above continued fractions are non-periodic (if there exist re- 
ally infinity of twins, Mersenne primes etc) they can not be solutions of a polynomial 
equations with rational coefficients of the degree 2, but we believe this statement 
remains true for rational polynomials of all degrees. Namely we are convinced that 
all considered above continued fractions u,U2,Ud,UFi,Uq,UM,Ur± are transcenden- 
tal, however we were not able to prove it and this problem seems to be extremely 
difficult. But if say u or U2 is not transcendental what the particular polynomial 
equation with very special (mysterious) integer coefficients should it satisfy? 

It is well known that the Champernowne constant [1A\ built by concatenating 
consecutive numbers in the base b is transcendental: 



where (7^)?, denotes number k expressed in the base b (e.g. in the common in 
computer science notation the twelfth number in the hexadecimal base 6 = 16 is 
denoted C). In the human base 6 = 10 the Cio is given by: 

Cio = 0.12345678910111213141516171819202122232425262728293031 . . . 

Transcendentality of Cb is the corollary from the theorem proved by Kurt Mahler 
in paper [33] published in 1937. In fact in this paper [33] Mahler has proved more 
general result: the number a obtained by concatenating the values of the positive, 
integer-valued increasing polynomial f{k) in the base b: 




(73) 



a = /(l),/(2)J(3),... 
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is transcendental, where f{k)h denotes the digits of the value of f{k) in the base 
system b. The case of Champernowne constant is not mentioned in [33] explicitly 
but it follows for f{k) = k. Let us remark, that the continued fraction expansion 
of Cio behaves very erratically, with sporadic partial quotients of enormous size, 
for example the 19-th term is of the order 10^®^, what is the typical behaviour for 
the Liouville numbers, i.e. such numbers that for each n there will be infinity of 
rationals A/B such that |Cio ~ A/B\ < l/S", see spikes in the Fig. 8. 

The number Cce obtained by concatenation of 0. with the base 10 representa- 
tions of the prime numbers in order 

Cce = 0.235711131719232931374143475153 . . . 

is known as Copeland-Erdos constant [15]. In this paper Arthur Herbert Copeland 
and Paul Erdos have shown that Cce is normal, but apparently it is not proved 
that Cce is transcendental. 

We have mentioned in the Sect. 3 that s = [0; 1, 2, 3, 4, 5, 6, . . .] is transcendental. 
Thus we have the correspondence s Cb and u ^ Cce, where both elements of the 
former pair are shown to be transcendental and both members of the latter pair are 
conjectured to be transcendental. Of course we have p„ > n, tn > n, qn > n etc. but 
we do not know how the transcendence of u, U2, Uq follows from these inequalities. 

One of the transcendence criterion is the Thue-Siegel-Roth Theorem, which we 
recall here in the following form: 

Thue-Siegel-Roth Theorem: If there exist such e > that for infinitely many 
fractions the inequality 



holds, then r is transcendental. 

Let us stress, that e here does not depend on n — it has to be the same for all 
fractions An/Bn- This theorem suggests the following definition of the measure of 
irrationality /i(r): For a given real number r let us consider the set A of all such 
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< 



n 



1,2,3, 



(74) 



exponents 6 that 



< 



P 

Q 



< 



(75) 



has at most finitely many solutions (P, Q) where P and Q > are integers. Then 



jj{r) = inf^gA ^ is called the irrationality measure of r (sometimes any 6 fulfilling (75) 
is called irrationality measure and then the smallest 5 = /i is called the irrationality 
exponent). If the set A is empty, then /i(r) is defined to be oo and r is called 
a Liouville number. If r is rational then /i(r) = 1 and if r is algebraic of degree 
> 2 then /i(r) = 2 by the Thue-Siegel-Roth Theorem. There exist real numbers of 
arbitrary irrationality measure 2 < ^ < oo. Namely, the number 



[aj + 



(76) 



where a > 1, b = fi — 1, has the irrationality measure /i, see [TT]. For the constant 
e = lim„_5.oo(l + 1/n)"', which has the continued fraction of a regular form: 



e = [2; 1, 2, 1, 1, 4, 1, 1, 6, 1, 1, 8, 1, 1, 10, 1, 1, 12, 1, 1, . . . , 1, 1, 2n, 



(77) 



it is known that fi{e) = 2, see [SI pp. 362-365]. For tt it is known that 2 < yu(7r) < 
7.6304, see [10], and it is conjectured [21 p. 203] that fi^n) = 2. There is a bound 
6{n) > 2 for infinitely many n following from the fact that of any two consecutive 
convergents to r at least one satisfies the inequality 



Pn 
Qn 



< 



2Qr 



(78) 



see [221 Theorem 18] or [TT, Theorem 183] and further: of any three consecutive 
convergents to r, one at least satisfies 



Pr 



< 



(79) 
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see Theorem 20] or [271 Theorem 195]. Thus writing for convergents satisfying 



(78) or (79) appropriately e(n) = log(2)/ logQn and e'{n) = log(5)/2 logQ„ the in- 
equahty appearing in the Thue-Siegel-Roth Theorem will be satisfied for a given spe- 
cific n. Of course fractions Pn/Qn constructed in this way will have lim„_>oo ^{n) = 0, 
because Qn increase monotonically and there will be no exponent of Qn on the r.h.s. 



of (74) strictly larger than 2 and common for all n. In fact, Khinchin [29] has proved 
that almost all reals r have /i(r) = 2. 

The partial quotients of u, Ud, Ug, upi grow too slow to use the Davenport — Roth 
Theorem, but if the behaviour of the Mersenne primes ~ 2^"" mentioned at 
the end of Sect. 6 is valid, then we obtain for large n 

~1 



Qn > 2^2* , c= — = 2.101893933 . . . (80) 



and transcendence of u_m will follow from the Davenport-Roth Theorem (46). We 



illustrate the inequality (80) in the Figure 3 — the values of labels on the y-axis 



give an idea of the order of Qn{-M.n)'- the largest for n = 47 is of the order Q47 = 

gl.9984...xl0« _ 2.32928 ... X 1086789810l 

Usually the number r in question (for example e, vr, C(3), etc.) is given by some 
definition not involving continued fractions, but here we have expressions of u, . . . um 
only by continued fractions and we can not calculate directly the differences |r — 
Pn/Qn\, like it is possible for example for Liouville transcendental numbers or for e. 
For this last mentioned earlier, the possibility of explicit calculation of the 

difference |e — Pn/Qn\ gives that /i(e) = 2 see [B, pp. 351-371]. We do not have any 
ideas now how to express qa, f, u, Ud, Uq, ... independently by means of formulas not 
involving continued fractions. Nevertheless we have made the plot of the exponent 
6{n) in the difference: 



Qn 



^ n= 1,2,3,..., (81) 



where U stands for u, Ud, Ug, ... and Pn/Qn are convergents of continued fractions 
for U — it is a well known fact that convergents of continued fractions are the best 
rational approximations. 
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Fig. 3 Illustration of the inequality (80) for 3 < n < 47. Although the last points seem to 
coincide in fact Q47 = 2.32928 . . . x io867898io^ ^i^j^g 2^2«--^ ^ 1.21513 . . . x 1082034318 _ 

hundreds thousands orders of difference! 



In the Figures 4-8 we present plots of 6{n) for u, Ug and for vr, e as well as for Ciq 
for comparison. First we have calculated u, Ug, ... with 150000 digits accuracy from 
the generic definition by constructing the continued fractions with a many thousands 
partial denominators. Next we have calculated P„ and Qn for n until the difference 
\U — Pn/Qn\ was zero in prescribed accuracy. From the differences \U — ^\ we 
calculated the values of 5{n) = — log \U — ^\/ \og{Qn) and the sample of results is 



plotted on Fig. 4 and 5 for u and Ug. The bound following from (78) is fulfilled for 



all n < 2100. In the next Figures we present the plot of 6{n) for vr (Fig.6), e (Fig. 7) 
and Cio (Fig. 8). For Cio we have plotted 6{n) — 2 because values of this difference 
changes by many orders, in contrast to smooth behavior seen in the Figs. 4-7. The 
spikes seen in the Fig. 8 are similar to the behavior of the Liouville transcendental 
numbers, but the last statement in [33] asserts that Ciq is not the Liouville number. 
In jl5] J. Sondow has proved that: 



/i(r) = 1 + limsup 



logQ 



n+l 



2 + lim sup 



log a 



n 

From this we have for u as a„ = p„ ~ nlog{n) and for large n Qr 



(82) 
ra" that 
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yu('u) = 2 and the same for U2, Ug, but for Mersenne primes we get from the Wagstaff 
conjecture: 

fx{uM) < 2 + 2^"' - 1 = 2.47477 . . . (83) 

But if there is only finite number of Mersenne primes (and hence finitely many even 
perfect numbers), then ij,{um) = 1- In the Fig. 9 we present the plot of (5(A^„;n) = 
— \og\uM — Pn/Qn\/ ^og{Qn) and indeed the values oscillate around 1 + 2^ ^ = 
2.47477.... 

We propose the conjecture which indicates the way to deduce the transcendence 
of some continued fractions from transcendence of another ones: 

Conjecture (t^t): Let a = [oq; ai, 02, . . .], lim„_>oo ctn = c>o, and (3 = [bo] &i, &2, • • 
where a„, 6„ G N. Suppose there exists such no that for all n > rio the inequahty 
bn > cL-n holds. If a is transcendental then /3 is also transcendental. 

The condition lim„_j.oo = 00 is necessary: if a„, is bounded, say a„ < A for 
Vn, A G N, then (3 = [A; 2A, 2A, . . .] = y/l + A"^. Also for transcendental bn the 
above conjecture probably is not true. When the Conjecture (*) will be proved it 
will suffice for our purposes to invoke the transcendence of the continued fraction 



s = [0; 1, 2, 3, 4, . . .] (49), as for all examples from Sect. 2 we have an > n (then also 
some examples from Sect. 3 will be transcendental, like / with a„ = n\ and F with 



5 Final remarks 

We have raised above some questions concerning the transcendence of continued 
fractions with partial quotients given by prime numbers of a few special forms. We 
hope that the experimental results reported above will stimulate further research in 
the field. 

Acknowledgement I would like to thank professors Boris Adamczewski, Jaroslav 
Hancl and Michael Waldshmidt for e-mail exchange. 



26 




Fig. 4 The plot of 5(n) (black) and the bound 2 + log(2)/ log (5„ (red) for u following from 



the (78) up to n = 22380 (computatuions were done in precision 150000 digits and the 
value of |n — -P22380/Q2238O I was zero with accuracy 150000 digits). For each n the bound 



( |78| ) (as weU as bound ([79])) is fulfilled. 
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Fig. 5 The plot of 5{n) (black) and the bound 2 + log(2)/log(5n for Uq following from the 



(78). For each n the bound (78) is fulfilled 
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Fig. 6 The plot of S{n) (black) and the bound 2 + log(2)/ log Qn (red) for vr following 



from the (78) up to n = 20997 (computatuions were done in precision 50000 digits 
and the value of [vr — -P20997/Q20997I was zero with prescribed accuracy). 
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Fig. 7 The plot of 6{n) (black) and the bound 2 + log(2)/log(5„ (red) for e following 



from the (78) up to n = 17365 (computatuions were done in precision 45500 digits 
and at this n the value of |e — P17365/Q17365I was zero). The periodic structure of 
the continued fraction expansion for e is clearly seen. 
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Fig.8 The plot of 5cio(^) ~ 2 (black) and the bound log(2)/logQ„ (red) for Cio 



following from the (78). Because of the weird behavior of the partial quotients of 



continued fraction expansion for Cio we have subtracted 2 from 5{n) and plotted the 
graph with the y axis in the logarithmic scale. After each extremely large partial 
quotient a„ there is an abrupt drop in the values of 5{n) and the bound log(2) / log Qn 
with accompanying spike for n — 1, see (|5]). It took almost 4 days CPU time to get 
data for this plot. Collecting data was done in a few separate runs with different 
precisions. Because the partial quotient 0525 > 10^^^^°° and Q527 > 10'^''^^^^ the 
calculations for 526 < n < 1708 was performed with 1,000,000 digits precision, see 
eq.([5]). We stopped aX n = 1707 because ciitos > 10"^^^^"^*^. Spikes of 5cio(^) many 
orders higher then neighboring values suggest that Cio niay be the transcendental 
number of Liouville type, but it in contradiction with the last statement of the paper 

ESI. 
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Fig.9 The plot of -log \um - Pn/Qn\/ log{Qn) (black) and the bound 2 + log(2)/logQ^ 



(red) following from the (78) for 3 < n < 43. Here the value of um was obtained from 



all 47 known Mersenne primes with more than 120 millions digits: the accuracy was 
better than IQ^^^^^^^^^"^. The denominators Qn grow very fast and the bound 
2 + log(2)/log(Q„) tends quickly to 2. It took 12 days CPU time on the AMD Opteron 
2700 MHz processor to collect data for n < 40: the point n = 40 needed precision of 
almost 40,000,000 digits, as \um - Ao/(34o| = 1-5033 x 10-38789567^ ^j^j^g 
l/Ql^ = 4.501 ... X io-3i553835_ Jo Calculate the difference \um - Pn/Qn\ for 
n = 41,42,43 the precision of 100000000 digits was needed and one point took 6 
days on the same processor, as for example \um — -P43/Q43I < IQ^^^'^™'^'^'^ . 
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